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ABSTRACT 

We introduce a jet shape observable defined for an ensemble of jets in 
terms of two-particle angular correlations and a resolution parameter R. 
This quantity is infrared and collinear safe and can be interpreted as a scal- 
ing exponent for the angular distribution of mass inside the jet. For small 
R it is close to the value 2 as a consequence of the approximately scale 
invariant QCD dynamics. For large R it is sensitive to non-perturbative 
effects. We describe the use of this correlation function for tests of QCD, 
for studying underlying event and pile-up effects, and for tuning Monte 
Carlo event generators. 
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1 Introduction 



Over the past years, experimental and theoretical advances have made it possible 
to ask increasingly detailed questions about jets. Recently, there has been consid- 
erable interest in jet substructure pp. Much work has focused on identifying jets 
initiated by boosted heavy objects such as the Higgs boson and the top quark. A 
variety of techniques have been proposed to identify and characterize substructure on 
a jet-by-jet basis. It is natural to draw upon these techniques to motivate interesting 
observables for the study of QCD. In this paper we explore an observable defined on 
ensembles of jets. 

A natural candidate for such an observable is the ensemble average of the angular 
structure function introduced in Ref. [2]. We will see that this ensemble average has a 
clear physical interpretation in terms of an average scaling exponent. Its leading order 
behavior can be found from a napkin-sized computation. Below we will argue that 
this ensemble average provides an interesting observable for at least three different 
reasons. First, it is an infrared and collinear safe observable that, in the perturbative 
regime, measures the extent to which QCD jet dynamics is scale invariant. Second, 
because it is formulated in terms of two-particle correlations, it asks a particularly 
detailed question about jet substructure. We find that different Monte Carlo event 
generators give significantly different predictions. Consequently measuring this en- 
semble average could give valuable feedback on the performance of the Monte Carlo, 
with any disagreements pointing towards the need for additional tuning or improve- 
ment of the physics modeling. Finally, the ensemble average has a simple dependence 
on uncorrelated radiation, such as might be expected from underlying event (UE) 
or pile-up (PU) contributions to jets. This suggests that measurements of the en- 
semble average could yield useful information about the average contribution of the 
underlying event and pile-up to hard perturbative jets. More broadly, contributions 
to the ensemble average with different scaling behaviors will be more or less impor- 
tant at smaller or larger angular scales. In this sense the ensemble average exhibits 
a characteristic sensitivity to both perturbative and non-perturbative physics, with 
the former dominating at small angular scales and the latter becoming important at 
large angular scales. 

The outline of this paper is as follows. In Sec. [2] we review the definition of the 
angular structure function and introduce its ensemble average. In Sec. [3] we compute 
the leading order behavior of the ensemble average in the collinear approximation 
and discuss expectations for corrections to the leading order result. In Sec. [4] we 
investigate the sensitivity of the ensemble average to the underlying event and pile-up, 
formulating a procedure for measuring the average density of uncorrelated radiation 
for a given ensemble of jets. In Sec. [5] we gain additional insight into the physics 
of the ensemble average by considering ensembles of soft radiation in the transverse 
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regions of the detector. In Sec. [6] we discuss our results and present our conclusions. 



2 Average angular structure function 

It has long been appreciated that jets have a fractal-like structure. This point of 
view emerges naturally from the description of the parton shower as a probabilistic 
Markov chain. The authors of Ref. [3] have computed the fractal dimension of a jet, 
while Ref. |1] advocates the use of the fractal phase space introduced in Ref. [5J as a 
useful diagnostic tool for complex events. 

Correlation functions provide a convenient language for studying fractal systems. 
In particular they can be used to define fractal dimensions through their limiting 
behavior at small scales. With this in mind, let us review the pair of correlation 
functions introduced in Ref. j2]. The first is the 'angular correlation function,' defined 
a£] 

g(R) = J^PTiPTjARleiR - ARij) (1) 

where the sum runs over all pairs of constituents of a given jet and 0(x) is the Heav- 
iside step function. Here pxi is the transverse momentum of constituent i, and ARij 
is the Euclidean distance between i and j in the pseudorapidity (77) and azimuthal 
angle (0) plane: Ai??- = (7^ — r/j) 2 + (0j — 4>j) 2 . Infrared and collinear safety and 
z-boost invariance fix this as the unique form for a two-particle angular correlation 
function defined on the constituents of a jet, with the only remaining freedom being 
that the exponent of AR^ in Eq. [I] is arbitrary so long as it is positive. Q{R) is the 
contribution to a jet's mass from constituents separated by an angular distance of R 
or less. It is worth emphasizing that R does not mark the distance with respect to 
any fixed center. 

In the context of fractals, a correlation function c(R) gives rise to a corresponding 
correlation dimension T> c defined as [6]: 

V c = lira l °f c{R) (2) 

R->0 logi? V ' 

There is of course an immediate obstacle to using Eq. [T] and Eq. [2] to define the 
correlation dimension of a jet: the finite resolution of the detector makes the small 
R limit inaccessible. In addition, the fractal-like structure of a jet does not continue 
down to arbitrarily small scales, since, for a jet with transverse momentum p?, the 
parton shower is cutoff at an angular scale -R m i n > ^-qcd/pt- A sensible alternative 



^In Ref. [5] Q{R) is normalized so that G(R) — > 1 at large R; however, for the purposes of this 
paper it is convenient to leave S{R) unnormalized with dimensions of mass squared. 
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to Eq. [2] is to instead define an 'angular structure function' AQ(R) via a logarithmic 
derivative: 

Aa = rflog£ = RdG_ = & 

dlogtf GdR ^pTiPTjARlBiR- ARij) y ' 

For a jet with a finite number of constituents, the 5-function in Eq. [3] results in a 
noisy function of R. A convenient way to obtain a smooth version of AQ(R) is to 
replace the 5-function by a gaussian with a fixed width dR: 

/ x ex P (-x 2 /2dR 2 ) 
\J2-KaR 

In order to maintain the leftmost equivalence in Eq. [3] the O-function must also be 
replaced by an error function GdR with the same width dR. 

The angular structure function AQ (R) encodes the scaling of the angular correla- 
tion function at a particular value of R. In particular if Q(R) ~ R 13 then AQ(R) = (3. 
In this sense AQ (R) recovers a scaling exponent analogous to the correlation dimen- 
sion in Eq. |2} On a jet-by-jet basis AQ(R) exhibits dramatic peaks at prominent 
angular scales corresponding to separations between hard substructure in the jet. 
This property of AQ(R) is exploited in Ref. [2] to construct an efficient top tagging 
algorithm. In order to clearly observe scaling exponents, however, we will need to 
average over large ensembles of jets, since the number of final state particles in a 
single jet is too few to clearly observe fractal structure. Such an ensemble average 
will be the subject of the rest of this paper. 

We use the angular correlation function as our basic object, defining its ensemble 

average as 

1 N 

fc=i 

where N is the size of the ensemble and Q{R)k is the angular correlation function of 
the fcth jet. From this average, we define the average angular structure function: 

WW) - WKJ)^) 



D EfcLi T,i^ j PTk,iPTk, j ^R 2 i j s d R(R - AR^) 

N 

E fc =i Y J i^ j PTKiPTkjAR 2 ij Q AK (R - AR^) 



(6) 



3 



where 5dR,(-R) and 0dR(-R) are the gaussian and error functions with width dR, respec- 
tively. Note that the ensemble average is not defined as an average over N angular 
structure functions AQ(R) k : 

1 - 

(Ag(R))^-^2AG(R) k (7) 

k=l 

We make this choice for at least two reasons. First, the definition in Eq. [6] lends itself 
more easily to analytical computation and makes possible its interpretation as an 
average scaling exponent. Second, on an event-by-event basis, the angular structure 
function is quite noisy. Consequently, the ensemble average in Eq. [7] is significantly 
noisier than that in Eq. |6j 

Throughout this paper we will set dR = 0.04. Although nonzero dR sculpts the 
ensemble averages somewhat, especially near R = 0, for dR = 0.04 the effect is small 
enough that we need not consider it when calculating (G(R)) analytically. 



3 Calculating the average 



As we will now show, a striking property of (AQ(R)) is that its leading order 
behavior can be understood from a simple computation. This is in contrast to, e.g., 
the integrated jet shape ^(R) [Z], which requires a detailed calculation even for its 
leading order behavior [8]. This section is organized as follows. In Sec. 3.1 we compute 



the leading order behavior of (AQ(R)} in the collinear approximation. While it would 
be rewarding to perform a NLO computation of (AQ(R)), in the subsequent sections 
we limit ourselves to exploring some of the features we expect to emerge from a 
more complete calculation. This task will be made easier thanks to the clear physical 
interpretation of (AQ(R)) as a scaling exponent. First, in Sec. 3.2 we discuss the 
qualitative effect of the running of the strong coupling. Second, in Sec. |3.3|we explore 



higher order effects with an emphasis on the expected difference between quark and 
gluon jets. Finally, in Sec. 3.4 we briefly touch upon whether (AQ(R)) could be 
amenable to factorization. 



3.1 Collinear approximation 

To begin we compute the average value of the angular correlation function (G(R)) 
in the collinear approximation. To first order in a s , (G(R)) can be computed from 

a r R o do 2 r 1 

(G(R)) * ^P 2 T J o J dz P(z)z(l - z)6 2 Q(R - 6) (8) 
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where Rq is the radius of the jet algorithm and P(z) is the appropriate Altarelli-Parisi 
splitting function. Notice that, as discussed at the end of Sec. |2j for the purposes of 
this section it is enough to set dR = 0, although dR > will be needed for any actual 
measurement. Since we are interested in the interior of the jet, in all of the following 
expressions we will assume that R < Rq- This prevents us from making predictions 
about edge effects in (Q(R)), but we do not expect the collinear approximation to be 
a good approximation at larger R anyway. From Eq. [8] we find: 

{|CV quark jets 

(9) 
Y^C A + ^n F T R gluon jets 

We thus have the leading order result that the angular correlation function for 
QCD jets goes like R 2 . Consequently, for both quark and gluon jets, we have that 
(AQ(R)) = 2. Note that this holds for any jet algorithm. 

Some interpretation of this result is in order. The average angular structure func- 
tion (AQ(R)) is a measure of how energy is distributed within a jet. A typical QCD 
jet has a hard core with the structure of emissions in and around the core controlled 
by the soft and collinear singularities. The fact that for QCD jets (AQ(R)) = 2 at 
leading order tells us that QCD has a collinear singularity of strength d<d 2 /<d 2 . By 
contrast, if the energy were distributed uniformly over the entire cone of the jet, we 
would expect (Q(R)) ~ R A and (AQ(R)) ~ 4, at least up to edge effects. 

3.2 Running coupling 

The qualitative effects of including a running coupling are straightforward to un- 
derstand. We proceed by evaluating the strong coupling at a characteristic energy 
scale Q 2 = p 2 ^8 2 g(z), where g(z) is a function of the momentum fraction z that is 
regular for < z < 1. Typically, Q is taken to be the jet mass (g(z) = z(l — z)) or 
transverse momentum of the emitted gluon (g(z) = 1 — z), but for our purposes what 
is most important is that Q 2 is proportional to 9 2 . This implies that at small angles 
the strong coupling increases resulting in more radiation in the small angle region of 
phase space. Effectively, this increases the strength of the collinear singularity with 
respect to the fixed coupling expectation. Thus, we expect that the effect of a running 
coupling is to lower the value of (AQ(R)) with respect to the fixed coupling result. 

A simple calculation confirms this picture. Now including a running coupling we 

(S(R)) P(Mi - z) a '^ gW W - 0) (io) 

To lowest order the running coupling is 



QCD 
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where «o — 27r//3o and {3q is the leading coefficient of the QCD beta function. In the 
limit that p F R 3> Aq CD , the precise forms of g(z) and P(z) are irrelevant, and we 
find that 

(G(R)) oc f . (12) 



with the result that 



log 



(Ag(R)) ~ 2 - 



PtR 
Aqcd 

1 



log 



PtR 
Aqcd 



(13) 



13 



As expected, including the running coupling decreases the average angular structure 
function. To first order in 1/logf-^^M, Eq. 
jets. This effect is not negligible. For examp 
Aqcd = 300 MeV, we have that 1/log (-^ 



is true for both quark and gluon 
Cfor p T = 200 GeV, R = 1.0, and 
0.15. 



3.3 Higher order effects 

We can get an idea of the nature of the higher order corrections to the 0(a s ) fixed 



coupling result by continuing the calculation of Sec. 3.1 to 0{a s ). The kinematic 



identifications of the three parton final state are illustrated in Fig. [TJ which depicts 
the longitudinal momentum fractions and angles associated to each splitting. We 
impose angular ordering so that Q\ > 8 2 . The expression for the 0(a 2 ) contribution 
to (G(R)) is then given by: 



(G(R))c 



Of, 
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Pt 
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+z(l-z)(l-y) 



2k 



2tt 



{y% + 9{ + 2y6 l 6 2 cos. 



xe(R- yjy 2 6l + 61 + 2y6 1 6 2 cos (I?) } 



(14) 



Using unpolarized splitting functions and setting cos = results in a simple analytic 
formula given by 



(6(B)), 



as P 2 T R 2 \lc F +° 



2tt 



2tt 



16 Gf+ 120 GaGf+ 120 



n F T R C F ^j ^1+log^ 
(15) 
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Figure 1: Splitting diagram for the 0(a 2 ) contribution to (Q{R)). 
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(16) 



for gluon jets. From these expressions we can calculate (AQ(R)) to 0(a s ), finding 

7 



<A0(i2)>< 



for quark jets and 
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for gluon jets, where we have set n F = 5 in evaluating the color factors. Notice that, 
since 



d 



dlogR 



log [R 2 (l + elog(i? 2 ))] = 2 + 2e + 0(c' 



(19) 



(AQ(R)) q / g remain flat in R to this order in a s . In particular Eq. 19 implies that the 
u 2 R 2 terms in Eq. 15 and Eq. 16 do not contribute to Eq. 17 and Eq. 18 at 0(a s ). 
These terms do contribute to the normalization of the angular correlation function, 
but we do not expect them to be correctly given by the collinear approximation. We 

9 -ft 2 

suspect, however, that the R log ^§ terms are robust as far as the scaling exponent 
is concerned. 
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<Ag(R)> <Aff(R)> 




(a) Pythia8 (b) Herwig++ 



Figure 2: Average angular structure functions for ensembles of jets with pt > 200 GeV 
and no underlying event or initial state radiation. Red curves denote quark jets and blue 
curves denote gluon jets. These are anti-kT jets with jet radius Rq = 1.0 as generated with 
Pythia8 (left) and Herwig++ (right). See Appendix [A| for more details about the Monte 
Carlo. 



The increase of (AQ(R)) with respect to the leading order result can be loosely 
interpreted as resulting from a relative increase in the amount of energy radiated 
away from the center of the jet. This increase is largest for gluon jets as a result of 
the large associated color factors. 



Putting together the results from Sec. 3.1, Sec. 3.2, and Sec. 3.3 our expectations 



for the form of (AQ(R)} for quark and gluon jets are as follows. Apart from edge 
effects, we expect both (AQ(R)) q and (AQ(R)) g to be approximately flat in R as 
in the leading order result. Furthermore, if the effects of a running coupling are 
dominant then we also expect that (AQ(R)) q / g < 2. In addition we expect that 
(AQ(R)) g < (AQ(R)) g . Note that, even in the case when (AQ(R)) q / g are separately 
flat, an ensemble average over an admixture of quark and glue jets will in general not 
be flat. 

In Fig. [2] and Fig. [3] we plot (AQ(R)) for both quark and gluon jets as obtained 
from Pythia8 and Herwig++ dijet events. The ensembles are composed of anti-kT 
jets with R = 1.0 and px > 200 GeV. Initial state radiation (ISR) and the underlying 
event have been turned off. The Monte Carlo is roughly in accord with our expec- 
tations. Not surprisingly, our expectations are more in line with the angular ordered 
shower in Herwig++, where the expected difference between quark jets and glue jets, 
i.e. (AQ(R)) q < (AQ(R)) g , is unambiguous. The differences between the predictions 
for (AQ(R)) given by Pythia8 and Herwig++ are striking (see Fig. [3]), with (AQ(R)) 
being substantially higher for Herwig++ over a large range in R. This differencd^J 



* We have also performed Monte Carlo calculations of (AQ(R)) for an ensemble of jets pro- 
duced in the process e + e~ — > qq. In this case the differences between Pythia8, Herwig++, and 



s 




<A£(R» 




(c) Jets 



Figure 3: Average angular structure functions for three different ensembles of jets with 
Pt > 200 GeV and no underlying event or initial state radiation. These are anti-kT jets 
with jet radius Rq = 1.0 as generated with Pythia8 (solid) and Herwig++ (dashed). See 
Appendix [A] for more details about the Monte Carlo. 



which persists when the underlying event and initial state radiation are included (see 
Fig. 4(c) ), should be measurable and motivates making the measurement. It would 
be interesting to have a detailed understanding of how this qualitatively different be- 
havior emerges from the two codes; doing so, however, lies outside the scope of this 
paper. 



3.4 Factorization 



Before moving on we would like to draw a connection to a recent analysis on 
the factorizability of jet substructure observables in soft-collinear effective theory 
(SCET) [10J. Walsh and Zuberi determined necessary conditions for jet substructure 
observables to be factorizable in the sense that any such observable can be computed 



Pythia8+Vincia (showering done with Vincia 1.0.26 [H]) are small compared to the differences in 



Fig. 3(a) 
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as a direct product of universal terms up to power corrections. A central result of 
their analysis is that in order for factorization to hold an observable must not demand 
that soft modes individually resolve collinear modes. Because the angular correlation 
function is a two-particle correlation function, one might worry that it mixes soft and 
collinear modes in a way that upsets factorization. However, to leading power soft 
modes contributing to (Q(R)) do not resolve individual collinear modes: 

J2ptcPtsAR 2 cs Q(R - AR CS ) = P tjJ2p tsAR s^ R ~ AR ^ ( 20 ) 

c,s s 

Here the sum runs over collinear (C) and soft (S) modes and J refers to the jet. 
This suggests that the ensemble average of the angular correlation function should 
be factorizable in SCET. 



4 Effect of uncorrelated radiation 

In the previous section we discussed the shape of (AQ(R)) as determined by the 
perturbative final state shower. At a hadron collider the colored initial state means 
that the dynamics of jets cannot be understood separately from the underlying event 
and initial state radiation. For convenience in the following we will often collectively 
refer to any radiation that is not associated with the hard, perturbative final state as 
the "underlying event." In particular what we have in mind is comparably soft radia- 
tion that is uncorrelated with the hard scatter and which is approximately uniformly 
distributed in pseudorapidity and azimuth with transverse momentum density Atje- 



4.1 Background 

We would like to ask how this UE affects (AQ(R)) for an ensemble of jets of a 
given p T . Since by assumption A UE <C p T , we can neglect correlations of order A^ E . 
Furthermore, since most of the energy of the final state shower is localized in a hard 
core at the center of the jet, we can write the (9(Aue) contribution to (Q(R)) as 

(£(#)>a ue = PtAve r d<P I R'dR'R' 2 = % r A UE # 4 (21) 
Jo Jo 1 

This ansatz neglects edge effects due to the finite size of the jet, but these are expected 
to be small away from R = R for the approximately circular^ anti-kT jets used 
throughout this paper. The range of R for which this ansatz is valid will be smaller 

§See, e.g., Fig. 7 in Ref. [IT] for 'typical' jet shapes as generated by different algorithms. 
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for other sequential jet algorithms (e.g. Cambridge/ Aachen and kr) that yield more 
irregularly shaped jets. Thus we have: 



(G(R)) = (G(R)) P ert + ^PtAveR' (22) 

Since the perturbative piece goes approximately like R 2 , at large R the UE contri- 
bution is increasingly important. In the absence of the perturbative piece we would 
have (AQ(R)) ~ 4. Thus the inclusion of UE has the effect of increasing (AQ(R)) 
from its perturbative value of near 2 at small R towards the value of 4 characteristic 
of uniform radiation at large R. This behavior is evident in Fig. |4j Going to the 
average angular structure function, we can rewrite Eq. [22] as: 

(A ^ ))pert " (G(R)) - IptAue^ (23) 
Provided that we have some sort of estimate of the perturbative piece (AG(R)) pert , 



the distinctive i?-dependence of Eq. 23 can be used to fit for Aue- In the remainder 



of this section we will investigate a procedure for doing so. 



4.2 Procedure 

One possible procedure for extracting Aue for a given ensemble of jets works from 
the assumption that (AQ(R)) pert is approximately flat in R. We have seen some 
tentative evidence for this hypothesis in Sec. [3j We can estimate (AG(R)) pcrt by 
measuring (AQ(R)) at small R. Specifically we assume that the perturbative piece is 
exactly flat and given as 

(Ag(R)) pcrt = (Ag(R miQ )) = C (24) 



where -R m i n <C 1. We choose -R m i n = 0.25. With this ansatz, we can invert Eq. [23] to 
solve for Aue as a function of R: 

A (R) - 2{G{R)) {AG{R)) ~ ° (25) 
Aue(jR) " n PT R* 4 - C (25) 

The flatness of Aue(-R) over a wide range of R can then justify the ansatz a posteriori, 
and Aue(-R) can be averaged over a range in R to obtain Aue- Summarizing, the 
procedure would be: 

1. Measure C = (A£/(i? m i n )) and take this as an ansatz for the perturbative con- 
tribution to (AG(R)) for all R. 

2. Construct Eq. [25] for a range of R with R > i? m i n . 
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<Ag(R)> (Ag(R)> 



2.5 2.5 




°"8.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



(c) Pythia8 vs. Hcrwig++ 

Figure 4: Average angular structure functions for ensembles of jets with pj- > 200 GeV. 
The black curves have no underlying event nor any ISR. In contrast to the previous section, 
here underlying event and ISR are turned on for the colored curves. On the LHS, the 
Pythia8 samples make use of tune 4C, with the red curve having twice as much UE activity 
as the blue curve. On the RHS, the purple curve corresponds to Herwig++ tune LHC-UE7- 
2. For comparison, the bottom figure overlays the Pythia8 and Herwig++ curves. These 
are anti-kT jets with jet radius i?o = 1-0- See Appendix [A] for more details about the Monte 
Carlo. 



12 



Aue(R) 



0.5 



0.6 



0.7 



Figure 5: Aue(-R) curves as extracted via the procedure detailed in Sec. 4.2 The solid 



curves correspond to the Pythia8 samples in Fig. [4j while the dashed curve corresponds to 
the Herwig++ sample. The dotted curve corresponds to the same Herwig++ sample as 
the dashed curve, the difference being that the power law of the underlying event ansatz is 
changed from R 4 to i? 3 ' 3 . The matching was done at i? m i n = 0.25. 



3. Then Aue, the average value of Atje(-R) over a range in R, gives a measure of 
the average density of UE in the ensemble of jets. 

An additional sanity check on the extracted scales could be to measure Aue as a 
function of the number of primary vertices npy. If the above procedure makes sense, 
Aue should exhibit a clear linear dependence on npy, since pile-up will contribute 
uncorrelated radiation to the ensemble of jets in an amount proportional to npy — 1. 



4.3 Results 



Curves Aue(-R) as extracted from Monte Carlo data via this procedure are plotted 
in Fig.[5j For Pythia8, Aue(-R) is approximately flat over a broad range in R, justifying 



the ansatz in Eq. 22 Note that the scale of the red curve is about twice that of the 
blue curve. This is as expected, since the red sample has twice as much UE activity 
as the blue sample (but the same amount of ISR). The dashed Herwig++ curve, by 
contrast, is not flat and has an unambiguous downward slope. Changing the UE 
power law ansatz from R A to R 33 results in a Aue(-R) curve that is much natter, 
though still at a similar scale. That is to say that for Herwig++ the underlying event 
is not uniformly distributed, instead being clustered somewhat around the center of 
the hard, perturbative jet. This points to the possibility of fitting for the exponent 



in the UE ansatz of Eq. 22 but we do not explore this possibility any further. 

The average UE densities Aue defined from the extracted curves Aue(-R) are listed 
in Table [T] These are defined by averaging Aue(-R) i n the range from R — 0.4 to 
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Monte Carlo Sample 


-A-Trans 


Aue 


Pythia8 Tune 4C 


3.2 


GeV 


2.8 ±0.3 


GeV 


Pythia8 Tune 4C 


4.6 


GeV 


5.0 ±0.7 


GeV 


Pythia8 Tune 4C" 


6.0 


GeV 


7.2 ± 1.1 


GeV 


Herwig++ Tune LHC-UE7-2 


3.3 


GeV 


3.7 ±1.0 


GeV 


Herwig++ Tune LHC-UE7-2 (R 3 3 ansatz) 


3.3 


GeV 


3.2 ±0.4 


GeV 



Table 1: Table of extracted UE densities. A Trans is the average pt density in the transverse 
region and Aue is extracted from (AQ(R)). The 'error' bars quoted for Aue are the max- 
imum difference of the function Aue(-R) from its average, which is computed in the range 
from R = 0.4 to R = 0.8. Tunes 4C' and 4C" differ from tune 4C in that they have twice 
and thrice as much UE activity, respectively. See Appendix [A] for more details about the 
Monte Carlo. 



R = 0.8, with the perturbative contribution to (AQ(R)) matched at -R m i n = 0.25. 
Also listed are the average transverse momentum densities Axrans in the transverse 
region (for the definition of which see Sec. [5]). For all the samples we see that Aue 
tracks A-rrans as expected, although it tends to be somewhat larger. For the three 
Pythia8 samples we see how, as we scale the amount of UE activity, A Trans and Aue 
rise linearly and in tandem. 

We can ask how sensitive the extraction of Aue is to deviations from the flat- 
ness assumption. The largest dependence on the perturbative ansatz C in Eq. 25 is 



through the term (AQ(R)} — C. A large overestimation (underestimation) of C gives 
a large underestimation (overestimation) of the average UE density. For example, 
if for Pythia8 Tune 4C we let C be given by the (AQ(R)) obtained from Pythia8 
in the absence of UE/ISR, we find an 0(100%) correction to Aue- in contrast the 
corresponding correction for Herwig±± is only O(20%). 

The validity of the procedure explored in this section depends on a number of as- 
sumptions, and an assessment of its usefulness would require a detailed experimental 
study with a particular emphasis on the systematic uncertainties involved. In addi- 
tion, a better theoretical understanding of the perturbative contribution to (AQ(R)) 
would be needed to determine the degree to which the flatness assumption is war- 
ranted. Nevertheless, as the luminosity (and eventually the center of mass energy) 
continues to rise at the LHC, the increasing importance of UE and PU motivates the 
study of additional experimental handles on the impact of UE and PU on hadronic 
final state reconstruction. While something like the "jet-area/median" method intro- 
duced in Ref. [12] is much more ambitious, since it provides a handle on fluctuations 
in UE and PU, the procedure proposed here might be a useful complement to this 
and other methods. 
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5 Angular correlations in the transverse region 

To gain more intuition for the scaling information encoded in the average angular 
structure function, it is interesting to consider a qualitatively different ensemble of 
jets. Unlike many traditional jet shape observables, (AQ(R)) is a meaningful observ- 
able for regions of the detector populated with soft, diffuse radiation. In the following 
we will focus on the transverse region of dijet events, as traditionally employed in un- 
derlying event studies. We will explore three different models for the underlying 
event: (i) the analytically tractable Feynman- Wilson gas; (ii) a toy Monte Carlo that 
describes uniformly distributed mini-jets; and (iii) full Pythia8 simulation. 

The transverse region is defined with respect to the direction of the hardest jet 
or hardest charged track (0 = 0), with its two halves extending between = ±7r/3 
and <fi = ±27r/3 and rj = — ^ max and rj = ^ max (see Fig. [6J [13]. We choose to define 
the transverse region with respect to the hardest jet and set r/ max = 2.0. The two 
transverse regions are further distinguished by their total activity. The region with 
the largest scalar sum p? is denoted as the trans-max region and the other is the 
trans-min region. Historically, underlying event studies have often focused on these 
two regions, since they are less sensitive to contamination from the hard perturbative 
scattering |14j . In the following ensemble averages will always be defined on charged 
tracks in the trans-min region, although other variations are possible. 

The most commonly used models of the underlying event are based on a picture 
of hadronic collisions in which multiple parton interactions (MPI) contribute to the 
soft activity in each event [15J. These models proceed by extrapolating the 2 — > 2 
QCD matrix element to semi-perturbative energies. Spectator partons that are not 
involved in the hard collision are allowed to scatter, resulting in an approximately 
uniform population of low pt jets, at least up to moderate pseudorapidities. 
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5.1 Fey nman- Wilson Gas 



Following the MPI hypothesis, we will make a number of simplifying assump- 
tions and see what sort of (AQ(R)) we would expect to measure in the trans-min 
region. To begin with consider a model for the underlying event along the lines of the 
Feynman- Wilson gas |16j . We assume that the number of underlying event particles 
is distributed with a probability distribution p(n) that depends only on the number 
of particles n. Typically this is taken to be Poissonian. Also, we will assume that the 
transverse momentum of the underlying event particles has a distribution px{n) that 
is uniform in r\ and 0. 

In this case the average angular correlation function can be calculated from 

oo 

(S(R)) = J2 E PAn)p T (n)ARlQ dR (R - AR^) (26) 

n=2 i<j<n 

For fixed n, the sum over the pairs of constituents breaks up into Q) identical terms 
so that we can write 



@( R V = E (^)p( n )PT(n)PT(n) j [d V d4>] ARl 2 e dR (R - AR 12 ) 



(27) 



where [drjdcj)] = dr]idr] 2 d(j)id(j) 2 . Since all the R dependence is in the integral, taking 
the logarithmic derivative yields the simple expression: 

WW) - dlQ ^ R - n j [dr]d(j)] A R 2 l2 e dR (R - AR 12 ) { ] 

That is, assuming that the underlying event is uniformly distributed, we can compute 
(AQ(R)) without knowing the pt or particle number distributions. The average 
angular structure function becomes purely geometric. 

For a rectangle of size Arj x A<f>, (AQ(R)) can be computed for dR = 0. For 
example, for R < min(A?7, A0) we have: 

(AG(R)\ - 2nAr]A<t) - 4(Ar? + A<P)R + 2R2 - 1 8 A ^ + A< ^ p | Q( R i) ( 2q] 
{AQ{R)) ~ |A^A0 -I(A V + A^R + p " 4 " 5^^A^ + ° {R ] (29) 

The corresponding (AQ(R)) over the full range of R with dR = 0.04, A0 = tt/3, and 
Ar] = 4.0 is plotted as the black curve in Fig. [7j Edge effects quickly bring (AQ(R)) 
down to zero for R > 1. At small but nonzero R, (AQ(R)) approaches the value of 4 
that is obtained in the infinite-plane limit Arj, A<p — > oo. As R — > the smoothing 



pulls (AQ(R)} down from near 4, as would be the case for dR = (see Eq. 29). 
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<A£(R)> 
4f 




Figure 7: (AQ(R)) in the trans-min region. The black curve corresponds to the Feynman- 
Wilson gas. The purple curve corresponds to Pythia8 dijet events with prjet > 200 GeV. 
The red curve corresponds to the toy Monte Carlo, described below. Finally, the blue curve 
corresponds to an ensemble of single "DLA mini-jets" as employed in the toy Monte Carlo. 
See Appendix [A] for more details about the Monte Carlo. 



5.2 Toy Monte Carlo 

With Monte Carlo we can explore the (AQ(R)) that results from underlying event 
models more complicated than the simplest Feynman- Wilson gas. In particular we 
would like to understand the (AQ(R)) that results from MPI models. We use the 
following toy Monte Carlo to test whether we understand the key physics ingredients 
entering into the form of (AQ(R)) obtained from Pythia8. First, we generate n mini- 
jets, where n is selected from a Poisson distribution with mean (n) = 0.30 per unit 
area of azimuth/pseudorapidity. We choose the transverse momenta of the mini-jets 
from the distribution 



P(PT) ~ - A 2 ) (30) 

Pt 

where A = 1 GeV. To generate the emissions that make up each mini-jet we employ 
a sort of double logarithmic approximation (DLA) scheme that mimics the parton 
shower. That is, for each mini-jet we generate a fixed number of emissions, N g = 10, 
where each "gluon" is generated according to the distribution 

/ „ , s d9 dz dd> . . 

where 9 is the angular separation between the gluon and the center of the mini-jet, z 
is the energy fraction of the gluon, and <ft is the azimuthal angle around the mini-jet 
axis. The 9 and z distributions are cut off at finite values, with 0.1 < z < 1.0 and 
0.01 < 9 < 1.2. 
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Figure 8: (AQ(R)) for ensembles excluding events where the pr of the hardest charged 
particle in the trans-min region exceeds a given py mEO! . (AQ(R)) for the Feynman- Wilson 
gas is shown in black, while the colored curves correspond to Pythia8 dijet events with 
PTjet > 200 GeV. The purple curve has no px cut, while the blue, green, orange and red 
curves correspond to prmux of 16 GeV, 8 GeV, 4 GeV and 1 GeV, respectively. See Appendix 
[A] for more details about the Monte Carlo. 

The (AQ(R)) for an ensemble of single "DLA mini-jets" is plotted as the blue curve 
in Fig. [7| By construction, for R < # max = 1.2 DLA mini-jets have an average angular 
structure function near the perturbative value of 2. Since MPI models are dominated 
by two scaling behaviors, (i) the perturbative R 2 scaling of the substructure of the 
mini-jets at small R and (ii) the Feynman- Wilson scaling at large R, we expect our 
toy model to yield a (AQ(R)) between the limiting (AQ(R)) curves of (i) and (ii). 
This is what we see in Fig. |7[ where the red curve interpolates between the blue 
and black curves, with a smooth transition between the two regimes at intermediate 
R. Furthermore, we see that the toy Monte Carlo does a good job of describing the 
average angular structure function obtained from Pythia8. 

5.3 Emergence of jets 

In Fig. [8] we illustrate the effect on (AQ(R)) when the ensemble is defined to 
exclude events where the pr of the hardest charged particle in the trans-min region 
exceeds a given pTmax- The Feynman- Wilson curve is shown for comparison in black. 
As the cut is decreased from infinity, (AQ(R)) approaches the Feynman- Wilson curve 
at small R. This is because the px cut has the effect of removing events with harder 
(and therefore jettier) structure in the transverse region. The soft events that remain 
are composed of soft mini-jets whose substructure is less jetty due to the smaller 
dynamical range of the parton shower. This effect should be measurable in data. 
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6 Discussion and conclusions 



The angular correlation function provides a direct probe into the scaling behavior 
of QCD. Because it is formulated in terms of two-particle correlations, it asks a 
particularly detailed question about jet substructure. Many jet shape observables 
(the integrated jet shape, angularities, etc.) are linear in the momenta of the jet 
constituents and are explicitly defined with respect to the jet center, i.e. they are radial 
moments of one kind or another. In this sense they can be thought of as accessing 
one-particle correlations, and the measurement of (AQ(R)) would be expected to 
provide orthogonal information about the substructure of jets. 

This comment is especially relevant to the observed differences in (AQ(R)) be- 
tween Herwig++ and Pythia8. If these event generators are being tuned to exper- 
imental data for observables like the integrated jet shape, then it is not surprising 
that they should give different predictions for two-particle correlations. Understand- 
ing these differences will be important if measurements of (AQ(R)) are to be used 
for improving Monte Carlo tunes. For example, the evolution variable of the par- 
ton shower affects both the structure of emissions in the jet as well as the scale at 
which the running coupling is evaluated. In addition, numerous other effects such as 
hadronization, color reconnections, and the details of the UE model will enter into 
the prediction for (AQ(R)). One would also like to understand how matched samples 
of jets, with their different treatment of hard, wide-angle radiation, affect the form 
of (AG(R)). 

An especially interesting property of (AQ(R)) is that it can be interpreted as an 
average scaling exponent. This makes the leading order result (AQ(R)) = 2 partic- 
ularly simple to understand and provides useful intuition for the higher order effects 
explored in Secj3] It also forms the basis of the two applications of (AQ(R)) ex- 
amined in Sec. [4[and Sec. [5] In the first case, reliable extraction of Aue will require 
a better theoretical understanding of the perturbative contributions to (AQ(R)). In 
particular it remains to determine to what degree the flatness assumption made in 
Sec. |4.2 is warranted. In the second case, forming the ensemble average in the trans- 
verse region provides additional insight into the physics encoded in (AQ(R)). We find 
that MPI models predict a quasi-universal form for (AQ(R)) with (AQ(R)) ~ 2 at 
small R and the large R form following the Feynman- Wilson gas. Although the jetty 
nature of UE is already well established, a measurement of (AQ(R)) in the transverse 
region has the nice property that it exhibits both the perturbative substructure and 
the uniform distribution of mini-jets. 

In this paper we have argued that the ensemble average of the angular structure 
function makes for a particularly interesting jet shape observable. From a theoretical 
point of view its interpretation as a scaling exponent is especially compelling. From 
an experimental point of view the possibility of measuring the average contribution of 
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the underlying event and pile-up to hard perturbative jets is intriguing. To go further 
will require experimental input, and we hope that (AQ(R)) might be measured at 
the LHC. 
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A Monte Carlo 

Here we give details of how the various Monte Carlo samples were generated. 
Throughout Herwig++ refers to Herwig++ v. 2.5.1 [T71 - I20] . and Pythia8 refers to 
Pythia8 v. 8.150 [22H2I]. To cluster jets, we use the Fast Jet v. 2.4.2 [25] implemen- 
tation of the anti-kT algorithm (26]. Note that the final state shower in Herwig++ 
is angular-ordered, while Pythia8 has a p^-ordered shower. We simulate pp collisions 
at a center of mass energy of Eqm = 7 TeV. No attempt is made to model detector 
effects, with particle-level information being used in all cases. In Sec. |3j the samples 
have been generated with MPI and ISR turned off. In Herwig++ this is done via the 
flags 

ShowerHandler :MPIHandler NULL 
SplittingGenerator : ISR No 

while in Pythia8 it is accomplished via the flags 
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PartonLevel : MI = off 
SpaceShower : QCDshower = off 

Also, in Sec. [3] we use ensembles composed of either quark or gluon jets. These are 
obtained by choosing the hardest jet in each event from samples generated with pure 
quark or glue final states, i.e. gg — > qq and qq — > qq in the case of quark jets and 
99 ~~ >* 99 an d qq — » gg in the case of gluon jets. For the event samples in Sec. [4] 
MPI and ISR are again turned on. For the underlying event Herwig++ makes use 
of tune 'LHC-UE7-2' [2]] and Pythia8 makes use of tune '4C [Ml EZj. In all cases 
we have checked that the pt distributions are similar between corresponding Pythia8 
and Herwig++ samples. Also, in Sec. [4] we use two Pythia8 event samples (tunes 4C 
and 4C") in which MPI activity has been increased by factors of 2.0 and 3.0. This is 
done by dialing the parameter Multiplelnteractions : Kf actor. 

Finally, in Sec. [5] the ensemble averages over the trans-min region come from a 
dijet sample where the px of the leading jet is greater than 200 GeV. The dijet sample 
is generated with Pythia8 using tune 4C, and the ensemble averages only make use 
of charged particles. Because of the large angular separations between charged tracks 
that can occur in the trans-min region, the ensemble average is modified to include 
a -Rmin prescription. That is, a given event only contributes to the numerator and 
denominator of (AQ(R)} in Eq. [6] for R > R mm , where -R m i n is the minimal angular 
separation between charged tracks in the trans-min region of that particular event. 
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